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Abstract. It has been recently demonstrated experimentally that graphene, or single-layer carbon, is a 
gapless semiconductor with massless Dirac energy spectrum. A finite conductivity per channel of order of 
e 2 /h in the limit of zero temperature and zero charge carrier density is one of the striking features of this 
system. Here we analyze this peculiarity based on the Kubo and Landauer formulas. The appearance of a 
finite conductivity without scattering is shown to be a characteristic property of Dirac chiral fermions in 
two dimensions. 
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We start with the Hamiltonian of a two-dimensional 
gapless semiconductor 
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Graphene, or single layer carbon pQ, demonstrates unique 
electronic properties. It has been shown recently (2113 that 
the charge carriers in graphene are massless Dirac fermions 
with effective "velocity of light" of order of 10 6 ms -1 . 
Graphene provides unexpected connections between con- 
densed matter physics and quantum held theory; in par- 
ticular, a new kind of quantum Hall effect observed in and the corresponding expression for the current operator 
graphene, that is, half-integer quantum Hall effect |U|3J Q5| 

can be considered as a consequence of the famous 
Atiyah-Singer index theorem The latter guarantees the 
existence of macroscopically large number of chiral states 
with zero energy in external magnetic field. 



Another amazing property of graphene is the finite 
minimal conductivity which is of the order of the conduc- 
tance quantum e 2 /h per valley per spin; it is important 
to stress that this is a "quantization" of the conductiv- 
ity rather than of the conductance 2\ This is not only 
very interesting conceptually but also important in light of 
potential applications of graphene for ballistic field-effect 
transistors pQ. Therefore the physical origin of the mini- 
mal conductivity is worth the special consideration which 
is a subject of this Note. 

Numerous considerations of the conductivity of a two- 
dimensional massless Dirac fermion gas do give this value 
of the minimal conductivity with the accuracy of some 
factor of order of unity pmi^lF?llTnilTTUT^lFT^j . It is really 
surprising that in this case there is a final conductivity for 
an ideal crystal, that is, without any scattering processes 

This fact is important since without complete under- 
standing of the ideal crystal case one can hardly hope to 
have a reliable answer for the realistic case with disorder 
and electron-electron interactions. Here we use the Lan- 
dauer formula 14 to clarify the physical meaning of this 
anomaly. 



J = «>X>^p = X>. ( 2 ) 
p p 

where v is the electron velocity, a — (<r x , cr y ) are Pauli ma- 
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trices, p is the momentum, and ^ 



pseudospinor electron operators. Here we omit spin and 
valley indices (so, keeping in mind applications to graphene, 
the results for the conductivity should be multiplied by 4 
due to two spin projections and two conical points per 
Brillouine zone). Straightforward calculations give for the 
time evolution of the electron operators 
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and for the current operator 

j(t) =h(t)+h(t)+il(t) 
Jo(*) = -£*p^p 
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where e p = vp/Ti is the particle frequency. The last term 
in Eq.QJ corresponds to the "Zitterbewegung", a phe- 
nomenon connected with the uncertainty of the position 
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of relativistic quantum particles due to the inevitable cre- 
ation of particle-antiparticle pairs at the position mea- 
surement |16II17| . Classical models for this phenomenon 
are discussed, e.g., in Ref.[T%| and references therein. 

In terms of condensed matter physics, the Zitterbe- 
wegung is nothing but a special kind of inter-band tran- 
sitions with creation of virtual electron-hole pairs. The 
unitary transformation generated by the operator U p = 
1/V2(1 + im p a), where m p = (cos0 p , — sm(f> p ) and (f> p 
is the polar angle of the vector p, diagonalizes the Hamil- 
tonian H p = diag (—vp, vp) and thus introduces electron 
and hole states; after this transformation the oscillating 
term in Eq.Q corresponds to the inter-band transitions, 
e.g. 



rtj*U v = ev 



i sin D e^ p 2i£pt 
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To calculate the conductivity a (u>) we will try first to 
use the Kubo formula which reads for two-dimensional 
isotropic case: 



oo 



°(") = ^f J d\(j(t-i\)j) 
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where (3 = T^ 1 is the inverse temperature, A is the sam- 
ple area. In the static limit uj = taking into account On- 
sager relations and analyticity of the correlators (j (z)j) 
for — j3 < Imz < one has [T§| 



oo 

a = U J dt{ i {t)i) 
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Usually, for ideal crystals, the current operator commutes 
with the Hamiltonian and thus j (t) does not depend on 
time. In that case, due to Eq.JfJJ) the frequency-dependent 
conductivity contains only the Drude peak 



ctd (uj) — lim Qz)-5 (uj) 
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Either the spectral weight of the Drude peak is finite and, 
thus, the static conductivity is infinite, or it is equal to 
zero. It is easy to check that for the system under consider- 
ation the spectral weight of the Drude peak is proportional 
to the modulus of the chemical potential (cf. Eq.(44) 
of Ref.[13|) and thus vanishes at zero doping (/x = Q). It is 
the Zitterbewegung, i.e. the oscillating term ji it) which 
is responsible for nontrivial behavior of the conductivity 
for zero temperature and zero chemical potential (that is, 
the gapless semiconductor case). A straightforward calcu- 
lation gives a formal result 



oo 

,-£/**■(«) 
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where one delta-function originates from the integration 
over t in Eq.Q and the second one - from the derivative 



of the Fermi distribution function appearing at the cal- 
culation of the average over product of Fermi-operators. 
Of course, the square of the delta function is not a well- 
defined object and thus Eq.® is meaningless before spec- 
ification of the way how one should regularize the delta- 
functions. After regularization the integral in Eq.© is fi- 
nite, but its value depends on the regularization proce- 
dure. It is not surprising therefore that two different ways 
of calculations in Ref.[S] led to two different answers. Our 
derivation, at least, clarifies the origin of these difficulties: 
it is the Zitterbewegung, or, physically, the impossibility 
to localize ultrarelativistic particles and to measure their 
coordinates. 

At finite frequency and finite chemical potential the 
Zitterbewegung contribution to the expression © coin- 
cides with the result for inter-band conductivity found in 
Ref-HSI- 

Despite this derivation cannot give us a correct numer- 
ical factor, it opens new way to qualitative understanding 
of more complicated situations. For example, the minimal 
conductivity of order of e 2 /h per channel has been ob- 
served experimentally also for the bilayer graphene PU] 
with the energy spectrum drastically different from that 
for the single-layer case. The bilayer graphene is a zero-gap 
semiconductor with parabolic touching of the electron and 
hole bands described by the single-particle Hamiltonian 

[2WT] 



H , (p x -i P ) /2m 

APx + Wy) /2m 







(here we ignore some complications due to large-scale hop- 
ping processes which are important for a very narrow 
range of the Fermi energies |23)- The Hamiltonian H10J1 
can be diagonalized by the unitary transformation U p with 
the replacement <p p — * 2<j) p . Thus, the current operator 
after the transformation takes the form J^J with the re- 
placement v — * p/m, e~ % ^ v — > e~ 2l ^ p . In contrast with 
the single-layer case, the density of electron states for the 
Hamiltonian Hl()[) is finite at zero energy but the square 
of the current is, vice versa, linear in energy. As a result, 
we have the same estimation ||{5J|, with the accuracy of 
additional factor 2. 

To circumvent the problem of ambiguity in the expres- 
sion for a in Eq. JHJ we now follow the alternative Landauer 
approach. Let us assume that our sample is a ring of length 
L y in y direction; we will use Landauer formula to calcu- 
late the conductance in x direction (see Fig. 1). There 
is still an uncertainty in the sense that the conductivity 
turns out to be dependent on the shape of the sample. 
To have a final transparency we should keep L x finite. On 
the other hand, periodic boundary conditions in y direc- 
tion are nonphysical and we have to choose L y as large as 
possible to weaken their effects. Thus, for two-dimensional 
situation one should choose L x <C L y . 

In the coordinate representation the Dirac equation at 
zero energy takes the form 



(K x + iK v )ifo = 
(K x - iK y )^ 2 = 



(11) 
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Fig. 1. Geometry of the sample. Thick arrow shows the direc- 
tion of current, ipt (solid line) and ipb (dashed line) are wave 
functions of the edge states localized near the top and the bot- 
tom of the sample, correspondingly. 



where K t = —i-^-. General solutions of these equations 
are just arbitrary analytical (or complex conjugated ana- 
lytical) functions: 



4>2 



ipi ( x + iy) , 

ip 2 (x - iy) ■ 



(12) 



Assuming that the conductance is equal to cfj^- one finds 
the contribution to the conductivity equal to e 2 /(irh). Ex- 
perimentally 0, it is close to e 2 /h, that is, roughly, three 
times larger than our estimation. The same result has been 
found earlier in Ref.jS] by one of the ways of derivation 
(the other one gives, instead, a factor tt/8). Note also that 
for the case of nanotubes (L x 3> L y ) one has a conduc- 
tance e 2 /h per channel, in accordance with known results 

The result a = e 2 /(irh) per valley per spin is found 
here for the case of ideal crystal. If one calculates in the 
simplest "bubble" approximation the conductivity in the 
presence of weakly scattering impurities and then put T = 
// = it leads to the same value |4"ll5ir7lllllll2| . However, 
one can hope that more transparent physical understand- 
ing of the origin of finite conductivity in ideal crystals 
which is provided by the Landauer formula will be useful 
to consider more complicated situations, such as the case 
of bilaver [2T)|. 

I am thankful to Andreas Ludwig, Andre Geim, and 
Kostya Novoselov for valuable discussions stimulating this 
work. 



Due to periodicity in y direction both wave functions should 
be proportional to exp(ik y y) where k y = 2nn/ L yi n — 
0, ±1, ±2, .... This means that the dependence on x is also 
fixed: the wave functions are proportional to exp (±2?:nx / L y 
They correspond to the states localized near the bottom 
and top of the sample (see Fig. 1). 

To use the Landauer formula, we should introduce 
boundary conditions at the sample edges [x = and 
x = L x ). To be specific, let us assume that the leads are 
made of doped graphene with the potential V < and 
the Fermi energy Ep — vkp = — V$. The wave functions 
in the leads are supposed to have the same y-dependence, 
that is, ipi^ {x, y) — -01,2 (%) exp (ik y y) . Thus, one can try 
the solution of the Dirac equation in the following form: 



Note added: After this work was basically finished ( cond-mat /051233' 
revised version) I have become aware of a relevant work 
by J. Tworzydlo, B. Trauzettel, M. Titov, A. Rycerz, and 
C. W. J. Beenakker (cond-mat/0603315) where a similar 
result for the transmission coefficient <|14[1 has been ob- 
tained, with a bit different choice of boundary conditions. 
)They have found also a sub-Poissonian shot noise in ideal 
graphene similar to that in disordered metals which gives 
a beautiful example of the importance of electron Zitter- 
bewegung. 



e lk * x + re lfcxX , 
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-01 (X) = 
02 (x) = 

where sin0 = k y /kF,k x 

tions of continuity of the wave functions, one can find the 
transmission coefficient 



be kyX , 



< x < L x (13) 
x > L x 



k 2 . From the condi- 



\t(kyW 



cosh (k y L x ) — sin 



(14) 



Further, one should assume that kpL x 3> 1 and put <fi ~ 
in Ea. H14|) . Thus, the trace of the transparency which is 
just the conductance (in units of e 2 /h) is 



TrT = 



E 



i 



cosh 2 (k y L x ) ttL 2 



(15) 
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